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The large-scale magnetic fields we observe today in galaxies and galaxy clusters could be the result 
of a pure quantum effect taking place during inflation, to wit, the creation of particles (photons) 
out from the vacuum in a curved spacetime. We show that, whenever the conformal invariance 
of electromagnetism is broken during inflation, the actual magnetic field spectrum, in the classical 
limit, is given by Bk ~ k 2 fanu, where 1 is the number of created photons with wavenumber 

k. In particular, a scale-invariant magnetic field of order of 10 12 G can emerge in the simplest 
model of cosmic magnetogenesis, the one where the inflaton is kinetically coupled to the photon. 
Moreover, and contrarily to the general belief, we show that such a model is free from the so- 
called strong-coupling and backreaction problems. This conclusion follows, indirectly, from the 
observation that post-inflationary electric currents, which in the literature are incorrectly supposed 
to freeze superhorizon magnetic fields after inflation, are indeed vanishing on superhorizon scales 
due to causality. 


PACS numbers: 98.80.-k,98.62.En 


Introduction. - Long time ago, Turner and Widrow [l| 
advanced the idea that the large-scale magnetic fields we 
observe today in galaxies and galaxy clusters (for a re¬ 
cent review on cosmic magnetic fields, see i) could be 
nothing else that relics from inflation. They also realized 
that, since Maxwell theory is conformally invariant in 
conformal-invariant spacetimes, classical magnetic fields 
(as those we observe today) cannot be produced in stan¬ 
dard electrodynamics. From this, the need arose to intro¬ 
duce nonstandard conformal-breaking terms in the elec¬ 
tromagnetic Lagrangian, such as nonminimal couplings 
between electromagnetism and gravity. Some years later, 
Ratra Q proposed a model in which the inflaton cf> (the 
scalar field responsible for inflation) is kinematically cou¬ 
pled to the photon via an exponential coupling of the 
form e^F 2 . This model has been successful for two rea¬ 
sons. First, it is able to produce scale-invariant magnetic 
fields whose intensity (~ 1CP 12 G) can directly explain 
cosmic magnetic fields. Secondly, the kinematic coupling, 
which was just postulated by Ratra, has been showed, by 
Martin and Yokoyama Q, to naturally emerge in string 
theory. However, soon after the framing of the Ratra 
model in a well-motivated, high-energy particle physics 
model, Demozzi, Mukhanov and Rubinstein [B| pointed 
out that inflation-produced magnetic fields able to ex¬ 
plain cosmic magnetization are either in an unmanage¬ 
able strong-coupling regime during inflation or strongly 
back-react on the dynamics of the inflaton. This sort 
of “no-go theorem” for inflationary magnetogenesis not 
only invalidated the predictions of the very elegant and 
simple model proposed by Ratra, but also the results of a 
plethora of other generating mechanisms put forward in 
the literature since the appearance of the seminal paper 
by Turner and Widrow. 

However, there is a flaw in the arguments of all the 
above cited works: it is assumed that post-inflationary 
electric currents froze the inflation-produced magnetic 


field on all scales, including superhorizon scales which are 
those of astrophysical interest for cosmic magnetic fields. 
But this would imply a violation of causality since, as first 
pointed out by Barrow and Tsagas @, post-inflationary 
currents are generated by microphysical processes during 
reheating and are, then, always vanishing on superhori¬ 
zon scales. 

The aim of this paper is to reconsider the generation of 
primordial magnetic fields at inflation in the kinetically 
coupled scenario by studying the creation of photons out 
from the vacuum in an expanding universe, a very pe¬ 
culiar phenomenon of non-conformal-invariant quantum 
field theories in curved spacetime [7) • In this framework, 
we will fix the “causality flaw” and, at the same time, 
we will show that there is no “no-go theorem” for infla¬ 
tionary magnetogenesis in the Ratra model, which indeed 
fully accounts for cosmic magnetic fields. 

The model. - It can be shown that the most general 
Lorentz-, C-, P-, CP-, CPT-, gauge-invariant, minimally 
coupled electromagnetic Lagrangian in four spacetime di¬ 
mensions, which involves classical, homogeneous, back¬ 
ground fields, has the simple form £ em = f(fa)CM • Here, 
Cm = —\F 2 is the standard (free) Maxwell Lagrangian, 
and the coupling function f(fa) incorporates all possi¬ 
ble interactions of the photon with scalar fields fa , which 
can be either composite or, as the inflaton cf> in the Ratra 
model, elementary. In order to have a positive-defined 
electromagnetic energy, and to avoid the strong-coupling 
problem [5j, f{fa) must be a positive-defined quantity 
greater than unity. Finally, in order to recover standard 
electromagnetism and not to spoil, thus, the predictions 
of the standard cosmological model, f(fa) must tend to 
one after reheating. 

We restrict our analysis to the case of a spatially flat 
universe described by the line element ds 2 = a 2 (dr/ 2 — 
dx 2 ), where r) is the conformal time and a(rj) is the ex¬ 
pansion parameter normalized to unity at the present 
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conformal time 770. Due to homogeneity and isotropy of 
both metric and background fields, the coupling f(<pi) 
turns to be a function only of the conformal time. 

Particle creation. - We now quantize our model. A 
general characteristic of quantum fields in curved space- 
time is that, since Poincare group is not in general a 
symmetry of the spacetime, the vacuum state cannot 
be uniquely defined as in the case of quantum theory 
in Minkowski spacetime. Nevertheless, we can define an 
“adiabatic” in-vacuum state |0, in) such that it is annihi¬ 
lated by the annihilation operators, 10, in) = 0, for all 
photon comoving wavenumbers k and polarization A, and 
normalized as ( 0 , in| 0 , in) = 1 . The in-vacuum state con¬ 
tains no in-particles in the “in-region”, 77 —> — oo, which 
will represent the temporal region when inflation starts. 
Such a vacuum can be constructed by solving the equa¬ 
tion of motion for the photon polarization states and then 
by fixing the constants of integrations appearing in the 
general solution by matching the latter with the corre¬ 
sponding adiabatic solution for 77 —> — oo 1 Accord¬ 
ingly, we can construct the in-Fock space based on | 0 , in) 
by repeatedly applying the creation operator aj) 1 )*' on the 
in-vacuum state. The electromagnetic vector potential in 
the in-Fock space reads 


a o?, x ) = 


A—1 * 


d 3 k £k,A (in) a (in) zkx 

( 2 tt )3 J 2 k k ’ A fe ’ A 


+ H.c., ( 1 ) 


where £k,A are the standard circular polarization vec¬ 
tors. The vector potential and its conjugate momentum 
tv = /A, as well as the annihilation and creation opera¬ 
tors, satisfy the usual commutation relations (hereafter, a 
dot indicates differentiation with respect to the conformal 
time). The photon wave function A^ n ^(r]) is a solution 


of the equation of motion Ak t x — (// f)Ak t x + k 2 Ak,x = 
0, and it matches the corresponding adiabatic solution 
in the in-region. In order to have a self-consistent 
quantization, it must satisfy the Wronskian condition 

4 (in) 4 (in)* 4 (in) 4 (in)* __ „ , , 

A k.X A k,X A k,X A k,X — ZIK /J- 

Consider, now, an “out-region” defined by r) —> +oo, 
which will represent any temporal region after reheating. 
An observer in this region may define an out-vacuum 
state |0,out) and construct from it the appropriate out- 
Fock space. He can then expand the electromagnetic field 
in terms of the out-annihilation and out-creation opera¬ 
tors, with the photon wave function AjA^ (77) matching 
the corresponding adiabatic solution at 77 —> +oo. 

It is clear from the above discussion that the two Fock 
spaces based on the two different choices of the vacuum 
state are both physically acceptable and must be then 
related. In particular, there will be a relation between 
the in- and out-modes Aj^ and Aj 43 ^, as well as a re¬ 
lation between the in- and out-creation and annihilation 
operators. In order to find these relations, let us re-scale 
the electromagnetic field as = y/f/ 2 k A^™’ out ^. 


The in- and out-ip modes satisfy the equation of mo¬ 
tion VA-,a = U k ipk, a, with t 4 (? 7 ) = -fc 2 + V 7 /VT- 

If and ip^\ are two solutions of this equation, 

the following inner product is conserved, {ip^\\ipk\) = 
~ We can then introduce the 

time-independent quantities ctk,x = (V’^'a IV , i°A t '**) anc ^ 
Pk ,A = — Accordingly, we can expand 

the in -ip mode in terms of the out -ip mode as ip^x = 


«© AT) 


Pk,\ipk > '\^*, where we used the fact that 
= = L Consequently, 


we have A^ = a*, aA^‘ ) + , which is the 

wanted relation between the electromagnetic in- and 
out modes. A relation of this type is know as Bogol- 
ubov transformation and the quantities ak,x and fik.x are 
called Bogolubov coefficients. They satisfy the relation 
|afc,A| 2 — \Pk,x \ 2 = 1, which can be easily derived from 
their defining equations. To find the relation between 
the in- and out-creation and annihilation operators, we 
insert the Bogolubov transformation in Eq. © and com¬ 
pare the result with the the expression of A defined in the 
out-Fock space. We find = ak,x — ftk x a -k a- 

From the above equation, it follows immediately that 
the two Fock spaces based on the two choices | 0 , in) and 
| 0 ,out) of the vacuum are generally different. In par¬ 
ticular, the in-vacuum state will contain out-particles so 
long as Pk,X ± 0 : n k ,A = ( 0 , in|| 0 , in) = \fi k ,x \ 2 , 

where is the number operator in 

the out-Fock space. This situation is understood as the 
creation of particles out from the vacuum by the chang¬ 
ing expansion parameter, and then by the changing grav¬ 
itational field. (Note that particles are created in pairs 
with opposite momenta, 7r_ k ,A = 7i k ,A-) However, such 
a particle creation is forbidden, for symmetry reasons, 
in conformal-invariant theories, such as standard electro¬ 
magnetism in a Friedmann-Robertson-Walker spacetime, 
a result known as “Parker theorem” Q • Accordingly, we 
must have no particle production if / is a constant. 

Actual magnetic field. - What we observe today (in the 
out-region) is the in-vacuum expectation value (VEV) of 
the squared magnetic field operator. The latter is 


(0,in|B 2 (77,x)|0,in) 




X=J° 


k ^ >fc ’ A 


in), 


where a 2 B = V x A is the magnetic field operator, 
'Pk’x = 5 'P^\-^'kx \ 2 is the magnetic power spectrum, 
and V/f{rf) = fc 4 /2-7r 2 a 4 the corresponding spectrum in 
the free Maxwell theory. Inserting the Bogolubov trans¬ 
formation in Vk’x , we find 

d:?(i) = +*£T%>+ATA), (2) 
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where 

C'Cfe+u.), 

p(out,0) _ 1 pM|^(°ut),2 
fc,A 2 ^ fc,A I ^ 

^” ix) = + 1 ) cos(n fc , A + 20 M ), 

with H fc , A = Arg(afc, A /3fc iA ) and 0fc, A (b) = ArgAjUfA 
The term Jdkk~ 4 (^1°^’^ + '^1°^ °^') f° rma lly 

equal to (n k , A ,n_k lA ,out|B 2 (? 7 ,x)|nk, A ,n_k, A ,out), 
where |nk.A, n_k.A, out) is the state containing ti^a out- 
particles with wavenumber k and helicity A, and n^k.A 
out-particles with wavenumber —k and same helicity 
A, so that (nk lA ,n_ k , A ,out|iV^° A t) |nk,A,«-k,A I out} = 

nk.A + n_k.A- In particular, the term J2\ Jdkk~ 1 V^ > J t ’ 0 ^ 
can be view as representing the out-vacuum mag¬ 
netic fluctuations, since it is formally equal to 
(0, out|B 2 (r;, x)|0, out). The quantity 
comes from the interference between the ak,\A^^ and 

Pk^f* wave functions and is, in general, different 
from zero. The physical origin of this term resides 
in the fact that the in-vacuum state is not generally 
proportional to the state |nk,A, n_k,A, out), in which case 
the interference term would be zero. 

In the out-region, where A[° ut ^ = e ~ lkv , we get, using 
Eq. ©, 


pM = p M 


- + n k + \/ Uk(n k + 1 ) cos(fifc - 2 kg) 


(3) 


evolution of the magnetic field is adiabatic. From these 
considerations, and defining the magnetic field on a scale 

1/k as Bk(if) = \JY1 a ^feA^’ we Sad, I n the limit n k ^ 1, 

Bfc(bo) - (kk 2 y/n^, (4) 

where <Y = ^/l + cos^flfc — 2kg\)/'K. Below, for the Ra- 
tra power-law coupling, we will find that f Ik approaches 
either the constant value 7r or — tt in the classical limit 
n k 1. Consequently, is an order-one constant. 

The result expressed by Eq. © is very different from 
the one we find in the standard literature of inflation- 
produced magnetic fields Q- In fact, the usual reasoning, 
which has been recently criticized by Tsagas Q, is the fol¬ 
lowing. One evolves superhorizon magnetic modes from 
the beginning up to the end of inflation. After reheating 
is completed, let us say at the time r ) e , such modes evolves 
adiabatically up to today, so that Bfc(bo) — a 2 {Ve)Bk{g e )- 
This result is incorrect since it reposes on the physically 
incorrect assumption that magnetic modes are frozen into 
the plasma even if they live on superhorizon scales. 

Ratra power-law coupling. - In the light of the above 
results, let us now reconsider the generation of magnetic 
fields at inflation in the well-known and well-studied case 
of a power-law coupling function which arises in the Ra¬ 
tra model. Let us take then f(g) = /*(b*/b) 2p if V < be, 
and f(g) = 1 otherwise, where be is the time when infla¬ 
tion ends. Here, /* > 1, p is a real parameter, p* = rj e if 
p < 0, and ? 7 * = bi if P > 0, where bi is the time when 
inflation begins. The Bunch-Davies-normalized ipk^ and 
V4° ut) wave functions are easily found. They are 


(hereafter, we drop the subscript A labeling photon po¬ 
larization states, if not strictly necessary). 

Let us concentrate on electromagnetic modes that be¬ 
haves classically in the out-region, namely modes with 
large occupation number, n k ^ 1- This is because only 
these modes may eventually explain the presence of cos¬ 
mic magnetic fields today. So far, we have neglected ex¬ 
ternal electric currents j M since they are produced only 
after inflation, during the process of reheating. Such 
currents, since generated by microphysical processes, are 
still vanishingon superhorizon scales after reheating due 
to causality | 6 |. Inside the horizon, instead, and if the 
electromagnetic field is classical, they can be written as 
J/i = (0, — acr c E), where a 2 E = —A is the electric field 
and a c is the conductivity of the primeval plasma. The 
latter is very high in the early Universe, so that post¬ 
inflationary electric currents wash out any electric field 
and freeze the inflation-produced magnetic field, which 
evolves then adiabatically, V^ oc a~ 4 , from the end of 
reheating till today (we are neglecting, here, any possible 
effect of magnetohydrodynamic turbulence [1]). Accord¬ 
ingly, Eq. © is valid up to the (positive) time bj., when 
a mode k crosses inside the horizon after inflation. Such 
a time is roughly defined by kg | ~ 1. For g > g±, the 


, (A 

n 


’(b) = 


d“*A) = 


A h "v( x '>' 

Jin) „(in) 

k, 1 ix I fc, 2 p —ix 
V2k V2k ’ 

(out) 

hv ?( x ) 

1 ~ix 


b < be, 
b > be, 


s/2 k 


(out) 

-JWKpix), b A be, 
b > be- 


Here, h Vp {x) = \fir/2 e z * (.* +v v) ^/xH^ p {x), H^\x) is 
the Hankel function of the first kind, v p = |i + p|, 

x = —A 77 , and are constant of integrations. The 

latter can be found by imposing the continuity of 
and V4° ut) , and their derivatives, at the time be- We 


find c (™) _ c (°ut) 

lmu c k, 1 — c fc, 1 


= (V’i m Vi° Ut) *}|r ) =r )e = OLk and 
= -4T ; = -{^k n) \^k Ut) )\v=Ve = Pk. The par- 


-fc, 2 


tide number is then 
1 

«k = T 


r i ( 1 v 



(*^e ) tli/p — 1 ( X e ) 


(5) 


where x e = —kg e , and it is zero for p = 0 (constant cou¬ 
pling function). Modes that are well inside the horizon 
at the end of inflation, x e 1 , are not efficiently pro¬ 
duced, rik ~ x ~ 4 (hereafter, we neglect numerical factor 














4 


in asymptotic expansions). The opposite is true for su¬ 
perhorizon modes (x e <C 1), «k ~ x e Up . Moreover, 
we find 0*, = 7rsgn(2t' p — 1) + 0(x e ) for x e <C 1. The 
most interesting case for inflationary magnetogenesis is 
when i/ p = 3/2 (p = —2 or p = 1) [Eq. © reduces to 
Uk = l/4xg]. This case, in fact, corresponds to a scale- 
invariant magnetic field B k in the out-region. Moreover, 
for a scale of inflation M of about 10 16 GeV (and as¬ 
suming, for simplicity, an instantaneous reheating), the 
intensity of such a field today, Bk(r]o) ~5x 10 -13 G, has 
the right value to directly explain cosmic magnetization 
(see, e.g., (I3|). 

Backreaction and renormalization. - The above analy¬ 
sis is valid only if the created (electro)magnetic field does 
not appreciably back-react on the dynamics of the Uni¬ 
verse. After inflation, electric fields inside the horizon 
are washed out by the high conductivity of the primeval 
plasma, while the sub-horizon inflation-produced mag¬ 
netic field evolves adiabatically. As it is easy to see, if 
such a held has to explain the observed cosmic magnetic 
fields, then its energy is always subdominant with respect 
to that of the Universe. During inflation, the in-VEV of 
the electromagnetic energy has to be much lower than 
the energy density of inflation, (0,in|p em (77, x)|0,in) -C 
Pi n f. Assuming, for the sake of simplicity, a de Sit¬ 
ter inflation, we have /?; n f = M 4 = 3H 2 / (8irG), where 
H M is the Hubble parameter during inflation and G 
the Newton constant. The electromagnetic energy den¬ 
sity is p em = /(E 2 + B 2 )/2. Using Eq. (HJ, we find 
(0,in|p em |0,in) = J^°dkk~ 1 p^ n \r]), where the electro¬ 
magnetic energy spectrum is 

Pt ] = pf[\KW\ 2 + ( 6 ) 


and “±” correspond to \p\ >1/2 and \p\ < 1/2, re¬ 
spectively. (For the scaling-invariant case, p = — 2 or 
p = 1, Eq. (J6J reduces to p^ = p^[l + (1 — §)x -2 + 
|(1 — p)x~ 4 }.) For large wavenumbers, x >> 1, p^ re¬ 
duces to the electromagnetic energy spectrum in the free 
Maxwell theory, p ^ = x 4 H 4 /2 tt 2 , which is ultraviolet di¬ 
vergent. This is a typical outcome of quantum theory 
in curved spacetime which can be cured by renormaliza¬ 
tion. For superhorizon modes, x -C 1, we have, instead, 
pj! n ' 1 ~ x 5 ~ 2i/ IpI H 4 . Different schemes of renormaliza¬ 
tion exist in the literature. However, whatever scheme 
is adopted, it must not change, on physical grounds, 
the expression of the energy spectrum in the classical 
limit. Classicalization is indeed realized for superhori¬ 
zon modes during inflation (see below). Therefore, we 
must have that = p £ n) for x -+ 0, where p ( k ^ hys 

is the renormalized (i.e. physical) energy spectrum, and 
^fcphys = 0 f° r x * +°° i n order to have a finite in- 
VEV of the electromagnetic energy density. The above 


behaviors for /4p hys are confirmed, for example, in the 
“adiabatic renormalization” scheme ;n]- We conclude 


that backreaction on inflation is completely negligible for 
\p\ < 2. This leaves the scale-invariant case free from the 
backreaction problem. 

Classicalization. - Let us now work in the Heisenberg 
representation. Following 12), 13], we introduce the time- 
dependent annihilation operator 


= J d3x ~^j :' [fc^(f?,x) + in(r?,x)]e jkx . 

The re-scaled electromagnetic field SH = \fj A and its 
conjugate momentum II = fF ^ (///)’!' (and, in turns, 
bk,\ and 6 k A ), satisfy the usual equal-time commutation 
relations. The evolution of the annihilation and creation 
operators follow from the Heisenberg equation of mo¬ 
tions, whose general solution can be represented through 
a Bogolubov transformation as 

bk,x{v) = pkb^\{-oo) + v k b\ x (-oo), 

b -k,x(v) = /4 6 -k,A(-°°) + v k b k.A(-oo). 

The Bogolubov coefficients pkip) an d u k( r i) satisfy the 
Bogolubov relation \p k \ 2 — \vk\ 2 = 1- We can re¬ 
late pk and v k to A^ n \ Since bk.\(—oo) = a k ‘ A and 
&l kA (— 00 ) = a^ kA , a straightforward calculation gives 

Pk - K = V?4 m) and Pk + u k = *V7 A^ n) /k. We 
can now use the standard parametrization of the Bo¬ 
golubov coefficients in term of the squeezing parame¬ 
ter rkirf) (which, without loss of generality, we can as¬ 
sume to be positive-defined), the squeezing angle <p k (ri), 
and the phase 0k{rj): p k = e~ l6k coshr^ and v k , = 
e i{e k + 2 ip k ) s i n ] lrfc- j n this case, using the expressions that 
relate pk and v k to , we find p\^ = p ^ cosh 2 ?-fe. 

We want now to determine when the inflation- 
produced electromagnetic field classicalizes. To this end, 
we re-consider the Bogolubov relation in physical units, 
\pk\ 2 — |^fe | 2 = ft, where H is the Planck constant. In the 
classical limit ft —> 0 , the above relation is satisfied only 
for r k —> + 00 . Therefore, large values of the squeezing 
parameter indicate that the system is in classical regime, 
and we conclude that electromagnetic fields during in¬ 
flation behave classically if their energy density is much 
greater than the energy density in the free Maxwell the¬ 
ory. For the power-law case discussed above, classical¬ 
ization is then realized, as anticipated, for superhorizon 
electromagnetic modes. 

Smooth coupling. In real models of inflationary mag¬ 
netogenesis, both the expansion parameter and the cou¬ 
pling f(r]) are smooth ( C °°) functions of the conformal 
time. Moreover, in the in- and out-regions (77 —> =Foo), 
the spacetime is slowly varying in standard Friedmann 
universes, in the sense that lim^ioo ^ = 0 for all 
integers n. If also the coupling function satisfies this last 
property (and this is a very reasonable assumption), then 
the in- and out-vacuum states are vacua of infinite adia¬ 
batic order [ 7 }. This implies that the number of produced 
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particle goes to zero exponentially (more precisely, faster 
than any power of k) in the limit k —> +oo Q . 

A numerical analysis 11] performed by using coupling 
functions which smoothly interpolate between f(rf) oc 
ij~ 2p during inflation and the asymptotic value f(rj) = 1 
in the out-region, confirms all the results obtained in the 
discontinuous power-law case in the long wavelength limit 
(this is because long wavelength modes are not strongly 
affected by variations of the coupling function on small 
timescales). Moreover, it clearly shows an exponential 
decay of n k for k +oo, instead of the fc -4 law derived 
in the discontinuous case. (Such kind of decays can be ex¬ 
plained analytically in semiclassical approximation Q) 
Tensor and scalar curvature perturbations. - Before 
concluding, we mention the fact that, in the Ratra power- 
law model and when backreaction on inflation is negligi¬ 
ble, the electromagnetic contributions to the primordial 
tensor curvature perturbations and to the “power spec¬ 
trum” , “bispectrum”, and “trispectrum” of primordial 
scalar curvature perturbations are fully compatiblejwith 
recent comic microwave background observations Q 
Conclusions. - We have pointed out that previ¬ 
ous calculations of the actual, inflation-produced mag¬ 
netic power spectrum, which rule out the Ratra inflaton- 
plroton kinematic coupling model, are incorrect. This is 
because post-inflationary currents, which are supposed 
to freeze superhorizon magnetic fields after inflation, are 
indeed vanishing on superhorizon scales due to causality 
arguments. The correct picture is, instead, that inflation¬ 
ary fields remain frozen into the primeval plasma only 
after they re-enter the horizon after reheating. Relating 
the magnetic spectrum after inflation to the occupation 
number of photons created out of the vacuum during in¬ 
flation, we have indicated the physical correct way to 


evolve magnetic fields from the beginning of inflation up 
to today, through their re-entering the horizon. We have 
found that the scale-invariant magnetic field produced 
in the Ratra model directly accounts for cosmic mag¬ 
netic fields if the scale of inflation is of order of 10 16 GeV. 
Moreover, such a field does not back-react on inflation 
and it does not suffer from strong-coupling problems. 

We thank A. Marrone for useful discussions. 


* Electronic address: leonardo.campanelli@ba.infn.it 

[1] M. S. Turner and L. M. Widrow, Phys. Rev. D 37, 2743 
(1988). 

[2] K. Subramanian, arXiv: 1504.02311 [astro-ph.CO]. 

[3] B. Ratra, Astrophys. J. 391, LI (1992). 

[4] J. Martin and J. Yokoyama, JCAP 0801, 025 (2008). 

[5] V. Demozzi, V. Mukhanov and H. Rubinstein, JCAP 
0908, 025 (2009). 

[6] J. D. Barrow and C. G. Tsagas, Phys. Rev. D 77, 107302 
(2008) [Phys. Rev. D 77, 109904 (2008)]. 

[7] N. D. Birrell and P. C. W. Davies, Quantum Fields in 
Curved Space (Cambridge University Press, New York, 
1982). 

[8] T. Kahniashvili, A. Brandenburg, L. Campanelli, B. Ra¬ 
tra and A. G. Tevzadze, Phys. Rev. D 86 , 103005 (2012); 
L. Campanelli, Phys. Rev. D 70, 083009 (2004); Phys. 
Rev. Lett. 98, 251302 (2007); Eur. Phys. J. C 74, 2690 
(2014). 

[9] C. G. Tsagas, arXiv: 1412.4806 [astro-ph.CO]. 

[10] L. Campanelli, Eur. Pliys. J. C 75, no. 6, 278 (2015); 
Phys. Rev. Lett. Ill, 061301 (2013). 

[11] L. Campanelli and A. Marrone, in preparation. 

[12] D. Polarski and A. A. Starobinsky, Class. Quant. Grav. 
13, 377 (1996). 

[13] M. Giovannini, Phys. Rev. D 61, 087306 (2000). 




